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Abstract. We find using Monte Carlo simulations of the spin-1/2 2D square 
lattice nearest neighbour quantum Heisenberg antiferromagnet that the high 
energy peak locations at (7r,0) and (7r/2,7r/2) differ by about 6%, (7r/2,7r/2) being 
the highest. This is a deviation from linear spin wave theory wrhich predicts equal 
magnon energies at these points. 

The simplest model describing quantum antiferromagnets is the nearest neighbour 
quantum Heisenberg model. Among the class of materials which to a good accuracy 
can be described by this model are the undoped high-temperature superconductors 
where the strongly interacting spins are located on a two-dimensional square lattice. 
Although simple to formulate, the Heisenberg model is not exactly solvable in 
dimensions greater than one, and approximations or numerical calculations are needed 
to compare the predictions of the Heisenberg model to experiments. 

For the 2D S = \/2 Heisenberg antiferromagnet on a square lattice, the time 
independent properties arc well understood |[ ^ , and measurements of for instance 
the correlation length ^, ||, ^, Q agree very well with the theoretical predictions. The 
situation concerning the dynamics is more unclear, and there is a need for definite 
predictions to which experiments can be compared. In particular, recent neutron 
scattering measurements on Cu(DCOO)2-4D20 and La2Cu04 |^ directly probe the 
magnon dispersion between the two points (7r/2,7r/2) and (7r,0) on the Brillouin zone 
boundary. These two materials, which are both considered to be physical realizations 
of the 2D S' = 1/2 Heisenberg antiferromagnet, show respectively a 6% decrease and 
a 13% increase in the magnon energy between (7r/2,7r/2) and (7r,0). These results 
are in contrast to the linear spin-wave approximation of the 2D Heisenberg model, 
which predicts equal magnon-energies at these points. While these deviations from 
linear spin- wave theory could be due to additional terms in the Hamiltonian describing 
each of the materials, it is also possible that there are corrections to linear spin-wave 
theory. In this Letter we aim at clarifying the predictions for the S* = 1/2 Heisenberg 
antiferromagnet on a square lattice at high energies, in particular at the special points 
(7r,0) and (7r/2,7r/2) in the Brillouin zone. 

The linear spin-wave approximation which is the zeroth term in an expansion in 
the parameter 1/ S gives the magnon dispersion 



LOk^US^l-^l (1) 

where 7^ = (cos kx -I- cos ky)/2. The wave-numbers are measured in units of the lattice 
spacing. Note that 7 is zero both at (7r,0) and (7r/2,7r/2). 
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Figure 1. Dynamic structure function S{q,u>) at q=(7r/2,7r/2) (dashed line) and 
q=(7r,0) (solid line). T = 0.5J and L = 32. 



The effect of the first order correction to linear spin- wave theory is to renormalize 
uniformly the magnon dispersion by a factor Z = 1.158. While there is a question as 
to what extent one can trust the spin- wave expansion for S = 1/2 it has been argued 
from Monte Carlo measurements that the only correction to linear spin-wave 
theory is such a uniform renormalization of the spectrum. Based on measurements 
at low energies this renormalization constant is found to he Z = 1.183, for 5=1/2 
[PI . This is similar to the ID system, where the exact solution gives a uniform 
renormalization Z = 7r/2. In 2D however, there are works which contradict the use of a 
single uniform renormalization constant. Employing the Dyson-Maleev representation 
of spin operators, Canali et al . [|2| found that the magnon energy at (7r/2,7r/2) is 
about 2% larger than the magnon energy at (7r,0) . Expanding around the Ising limit, 
Singh and Gelfand JlSf found a shallow minimum in the dispersion around (7r,0) giving 
the magnon at (7r,0) about 7% less energy than at (7r/2,7r/2). In addition, an approach 
to the Heisenberg Hamiltonian starting from the 7r-flux phase, a state with short-range 
antiferromagnetic order, predicted a deep local minimum around (7r,0) in the Brillouin 



To clarify this issue we have calculated the dynamic structure function S{q,uj) 
using the quantum Monte Carlo loop algorithm ||l^ , which among other useful features 
operates in continuous imaginary time [|l6j. While the simulations are performed in 
continuous imaginary time, the measurements of the spin-spin correlation function 
were written to an array with typically 200 points in the imaginary time direction. 
Typically 10^ configurations were used both for equilibration and measurements using 
a single-cluster implementation of the loop algorithm, and all data points are averages 
of at least five independent runs. The focus on high energy peaks makes it sufficient 
to do measurements at intermediate temperatures, T J, thus avoiding the low 
temperature region where the loop algorithm performs poorly. 

To get real-time dynamics from the imaginary time data we employed the 
maximum entropy method |]l7| , with a flat prior. As we restrict ourselves to only 
calculating relative magnon energies at two points in the Brillouin-zone, this choice 
of a priori information should not be crucial, although we expect that the flat prior 
overestimates the peak widths. As a check on the continuation procedure we evaluated 
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Figure 2. Magnon energies for tlie two points (7r,0) (filled circles) and 
(7r/2,7r/2) (open circles) in the Brillouin zone at different temperatures. The inset 
shows the relative difference r = (<.L)(-,r,o) ~ '^(7r/2,ir/2))/'^(ir/2.7r/2) iii percent. 



the sum rules corresponding to the -1., 0. and 1. moments of the dynamic structure 
function |p^ . They were all within the error bars of the quantities, S{q,iuJn — 0)/2, 
S{q,T = 0) and — e(2 — cos(q2;) — cos(gj^))/3, respectively, which were extracted directly 
from the imaginary time data (e is the energy per site). A typical picture of the 
dynamic structure factor thus obtained is shown in Figure ^ 

Following the approach of Makivic and Jarrell jTot , we determined the magnon 
energy from the normalized first moment tUq of the relaxation function 



F{q,iu) = 2(1 - e-^^)S{q,uj) {(3oJx{q)y 



(2) 



Because there is nothing that breaks the spin rotational symmetry in our Monte Carlo 
calculation this relaxation function is an average over the transverse and longitudinal 
relaxation function. This average resembles closely what is measured in neutron 
scattering experiments. 

Computation of Uq for our smallest systems, 4x4, gives no significant difference 
between W(7r,o) ^-nd 1^(77/2,77/2) • This is in agreement with the exact diagonalization 
study on small systems of Chen et al . However, lattices with 8x8 sites show a 
clear difference between 1^(77,0) B^nd W(7r/2,7r/2) for all temperatures studied. Performing 
a finite size analysis for L x L-systems, L = (4,8,16,32), we find that our results 
are consistent with the finite size behaviour ]lO| LUq^ ~ t^g.oo + Aq/ , where Aq is 
weakly temperature dependent and of order 10, and ^(77/2.71/2) < ^(7r o)- The resulting 
magnon energies for the infinite size system are plotted in Figure ^ as functions of 
temperature. 

It is clearly seen that the magnon energy at (7r,0) is lower than the magnon energy 
at (7r/2,7r/2). Extrapolating the relative difference to zero temperature we find that 
the magnon energy at (7r,0) is about 6% lower than at (7r/2,7r/2). This is in rough 
agreement with the result obtained by expanding around the Ising limit [p^ . 

Comparing this result with the measurements, it is seen that the zone boundary 
dispersion in Cu(DCOO)2-4D20 is indeed accounted for by the correction to spin- 
wave theory. Evidence of the zone boundary dispersion has also been observed in 
Sr2Cu304Cl2 which host two interpenetrating square lattices of = 1/2 spins |po| . 



Letter to the Editor 



4 



La2Cu04 on the other hand show a significant deviation from our resuit (19%). This 
is most likely due to higher-order spin couplings in the Hamiltonian describing this 
system. 

A qualitative explanation for the zone boundary dispersion can be achieved by 
considering the Hubbard model which for large-[//i at half-filling is equivalent to 
the Heisenberg model. A reasonable ground state ansatz for the Hubbard model 
which has a low energy and respects time-reversal invariance is the tt-Aux state in 
which the electrons behave as if they were subjected to a magnetic field of flux tt per 
plaquette. In this state the Fermi "surface" is located at the points (7r/2,±7r/2) in 
the Brillouin zone, and as the magnons are particle-hole excitations their dispersion 
will have minima at (0,7r) and (7r,0) as well as at (0,0) and (tTjTt). This qualitatively 
explains the zone boundary dispersion, but implies at the quantitative level gapless 
magnons at (7r,0) and (0,7r). However, the 7r-flux phase has no antiferromagnetic long- 
range order, and so it is necessary to do a more refined ground state ansatz to get a 
quantitative explanation. By considering a ground state consisting of both a tt-Aux 
state and a spin density wave (SDW) state which does have a staggered magnetic 
moment, Hsu[0 found using a Gutzwiller projection technique that the total energy 
of this combined state could be decreased from the pure 7r-flux state. Using the random 
phase approximation he found a shallower zone boundary minimum than for the pure 
TT-flux state. Though still much deeper than what is observed in the present study, 
the depth of the zone boundary dispersion at (7r,0) depends crucially on how much of 
the ground state is an SDW, as a pure SDW state gives a dispersion identical to linear 
spin wave theory[pT[. It is at present unclear whether a more refined analysis would 
lead to a better quantitative agreement with our result, but we feel that attributing 
the zone boundary dispersion to properties of the 7r-flux state is a plausible, simple 
and attractive scenario. 

Until now, there has been an unsatisfactory situation, where different 
approximative analytical approaches disagreed on the existence of a zone boundary 
dispersion in the 2D S" = 1/2 Heisenberg antiferromagnet on a square lattice. Recent 
experiments on physical realizations of the model system have strengthened the need 
to resolve this question. By performing finite temperature quantum Monte Carlo 
simulations, we have established that indeed there is a zone boundary dispersion. In 
a system, which is otherwise well described by a uniform renormalization of linear 
spin- wave theory, the zone boundary dispersion is a remarkable quantum effect. 
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